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October 12th, 2006 

Abstract A significant new development in the study 
of Anomalous X-ray Pulsars (AXPs) has been the recent 
discovery by INTEGRAL and RXTE of flat, hard X-ray 
components in three AXPs. These non-thermal spectral 
components differ dramatically from the steeper quasi- 
power-law tails seen in the classic X-ray band in these 
sources. A prime candidate mechanism for generating 
this new component is resonant, magnetic Compton up- 
scattering. This process is very efficient in the strong 
magnetic fields present in AXPs. Here an introductory 
exploration of an inner magnetospheric model for up- 
scattering of surface thermal X-rays in AXPs is offered, 
preparing the way for an investigation of whether such 
resonant upscattering can explain the 20-150 keV spec- 
tra seen by INTEGRAL. Characteristically flat emission 
spectra produced by non-thermal electrons injected in 
the emission region are computed using collision inte- 
grals. A relativistic QED scattering cross section is em- 
ployed so that Klein-Nishina reductions are influential in 
determining the photon spectra and fluxes. Spectral re- 
sults depend strongly on the magnetospheric locale of the 
scattering and the observer's orientation, which couple 
directly to the angular distributions of photons sampled. 
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1 Introduction 

Over the last decade, there has been a profound growth 
in evidence for a new class of isolated neutron stars with 
ultra-strong magnetic fields, so-called magnetars that 
include Soft-Gamma Repeaters (SGRs) and Anomalous 
X-ray Pulsars (AXPs) . Such a class was first postulated 
as a model for SGRs by Duncan & Thompson (1992), and 
later for AXPs (Thompson & Duncan 1996). The AXPs, 
are a group of six or seven pulsating X-ray sources with 
periods around 6-12 seconds. They are bright, possess- 
ing luminosities Lx ~ 1 35 erg s" 1 , show no sign of any 
companion, are steadily spinning down and have ages 
t <> 10 5 years (e.g. Vasisht & Gotthelf 1997). The steady 
X-ray emission has been clearly observed in a number of 
AXPs (e.g. see Tiengo et al. 2002, for XMM observa- 
tions of IE 1048.1-5937; Juett et al. 2002 and Patel et 
al. 2003, for the Chandra spectrum of 4U 0142+61), and 
also SGRs (see Kulkarni et al. 2003 for Chandra observa- 
tions of the LMC repeater, SGR 0526-66). A nice sum- 
mary of spectral fitting of ASCA X-ray data from both 
varieties of magnetars is given in Perna et al. (2001). 
This emission displays both thermal contributions, which 
have kT ~ 0.5 — 1 keV and so are generally hotter than 
those in isolated pulsars, and also non-thermal compo- 
nents with steep spectra that can be fit by power-laws 
dn/dE oc E~ s of index in the range s = 2 — 3.5 . 

Flux variability in AXPs is generally small, suggest- 
ing that even the non-thermal components experience a 
moderating influence of the stellar surface, rather than 
some more dynamic dissipation in the larger magneto- 
sphere. Yet the recent observation (Gavriil & Kaspi 2004) 
of long-lived pulsed flux flares on the timescale of sev- 
eral months in AXP IE 1048.1-5937 resembles earlier 
reports (Baykal & Swank 1996; Oosterbroek et al. 1998) 
of modest flux instability. There are also correlated long 
term variations in X-ray flux and non-thermal spectral 
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index in the source 1RXS J170849.0-400910, as identified 
by Rea et al. (2005). Moreover, Kaspi et al. (2003) and 
Gavriil, Kaspi & Woods (2002, 2004) reported bursting 
activity in the AXPs IE 2259+586 and IE 1048.1-5937, 
suggesting that anomalous X-ray pulsars are indeed very 
similar to SGRs, a "unification paradigm" that is cur- 
rently gathering support, but remains to be established. 

The recent detection by INTEGRAL and RXTE of 
hard, non-thermal pulsed tails in three AXPs has pro- 
vided an exciting new twist to the AXP phenomenon. In 
all of these, the differential spectra above 20 keV are ex- 
tremely flat: IE 1841-045 (Kuiper, Hcrmsen & Mendcz 
2004) has a power-law energy index of s = 0.94 between 
around 20 keV and 150 keV, 4U 0142+61 displays an in- 
dex of s = 0.2 in the 20 keV - 50 keV band, with a steep- 
ening at higher energies implied by the total DC+pulscd 
spectrum (Kuiper et al. 2006), and RXS J1708-4009 has 
s = 0.88 between 20 keV and 150 keV (Kuiper et al. 
2006); these spectra are all much flatter than the non- 
thermal spectra in the < 10 keV band. Also, no clear 
tail has been seen in IE 2259+586, yet there is a sug- 
gestion of a turn-up in its spectrum in the interval 10- 
20 keV (Kuiper et al. 2006). The identification of these 
hard tails was enabled by the IBIS imager on INTE- 
GRAL and secured by a review of archival RXTE PC A 
and HEXTE data. These tails do not continue much be- 
yond the IBIS energy window, since there are strongly 
constraining upper bounds from Comptel observations 
of these sources that necessitate a break and steepen- 
ing somewhere in the 150-750 keV band (see Figures 4, 
7 and 10 of Kuiper et al. 2006). Interestingly, Molkov 
et al. (2004) and Mereghetti et al. (2004) also reported 
evidence for hard tails in SGR 1806-20, so that the con- 
siderations here are germane also to SGRs in quiescence. 

Explaining the generation of these hard tails forms 
the motivation for this paper, which presents an ini- 
tial exploration of the production of non-thermal X-rays 
by inverse Compton heating of soft, atmospheric ther- 
mal photons by relativistic electrons. The electrons are 
presumed to be accelerated either along open or closed 
field lines, perhaps by elcctrodynamic potentials, or large 
scale currents associated with twists in the magnetic field 
structure (e.g. see Thompson & Beloborodov 2005). In 
order to power the AXP emission, they must be pro- 
duced with highly super-Goldreich-Julian densities. In 
the strong fields of the inner magnetospheres (i.e. within 
10 stellar radii) of AXPs, the inverse Compton scatter- 
ing is predominantly resonant at the cyclotron frequency, 
with an effective cross section above the classical Thom- 
son value. Hence, proximate to the neutron star surface, 
in regions bathed intensely by the surface soft X-rays, 
this process can be extremely efficient for an array of 
magnetic colatitudes. Here, an investigation of the gen- 
eral character of emission spectra is presented, using col- 
lision integral analyses that will set the scene for future 
explorations using Monte Carlo simulations. This sce- 
nario forms an alternative to recent proposals (Thomp- 



son & Beloborodov 2005; Heyl & Hernquist 2005) that 
the new components are of synchrotron or bremsstrahlung 
origin, and at higher altitudes than considered here. The 
efficiency of the resonant Compton process suggests it 
will dominate these other mechanisms if the site of elec- 
tron acceleration is sufficiently near the stellar surface. 
This prospect motivates the investigation of resonant in- 
verse Compton models. 



2 The Compton Resonasphere 

2.1 Energetics 

The scattering scenario for AXP hard X-ray tail forma- 
tion investigated here assumes that the seed energization 
of electrons arises within a few stellar radii of the magne- 
tar surface. This can in principal occur on either open or 
closed field lines, so both possibilities will be entertained. 
The key requirement is the presence of ultra-relativistic 
electrons moving along B, with an abundance satisfying 
the energetics of AXPs implied by their intense X-ray 
luminosities, Lx ^ 10 35 erg/sec above 10 keV (Kuiper 
et al. 2006). The hard X-ray tails have luminosities that 
are 2-3 orders of magnitude greater than the classical 
spin-down luminosity E SD ~ 8n 4 B%R 6 /(3P 4 c 3 ) due to 
magnetic dipole radiation torques. Here B p is the sur- 
face polar field strength. This signature indicates that 
other dissipation mechanisms, such as structural rear- 
rangements of crustal magnetic fields, power the AXP 
emission (e.g. Thompson & Duncan 1995; 1996). Let n e 
be the number density of such electrons, (j e ) be their 
mean Lorentz factor, and e ra( j be the efficiency of them 
radiating during their traversal of the magnetosphere (ei- 
ther along open or closed field lines). Then one requires 
that Lx ~ e ra d(7e)?n e c 2 (47me-Rc c ) if the emission col- 
umn has a base that is a spherical cap of radius R c . 
This yields number densities n e ~ 3xl0 17 I/x,35Arad(7e) 
cm~ 3 for scaled luminosities £x,35 = Lx/10 35 erg/sec, 
if R c ~ 10 6 cm. Therefore large densities are needed, 
though not impossible ones, since optically thin condi- 
tions for the surface thermal X-rays prevail provided that 
(je) 3> 1 , and e ra d is not miniscule. 

Comparing en e to the classic Goldreich- Julian (1969) 
density p GJ = V.E/47T = — fi.B/(27rc) for force-free, 
magnetohydrodynamic rotators, one arrives at the ratio 

en e _ 4, 670 L X ^P ,^ 

for AXP pulse periods P in units of seconds, polar mag- 
netic fields £?i 5 in units of 10 15 Gauss, and cap radii 
i?6 m units of 10 6 cm. Large electron Lorentz factors 
of 7 e 3> 10 2 — 10 3 and their efficient resonant Comp- 
ton cooling (i.e. e rac i ~ 0.01 — 1 ) are readily attained 
in isolated pulsars with B ~ 0.1 (e.g. see Sturner 1995; 
Harding & Muslimov 1998; Dyks & Rudak 2000), and 
such conditions are expected to persist into the mag- 
netar field regime. For e ra d(7e) *S 10 3 and Re ~ 1 , 
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the requisite density n e is super-Goldreich- Julian, but 
not dramatically so. This situation corresponds, how- 
ever, to acceleration zone cap radii R c that are con- 
siderably larger than standard polar cap radii R9 ca , p for 
AXPs. Here 9 cap — arcsin{(27r R/Pc) 1 / 2 } for a pulsar 
of radius R and period P . For P = 10 sec, this yields 
-R^cap ~ 4.6 x 10 3 cm. Concentrating the relativistic elec- 
trons in such a narrow column yields charge densities 
far exceeding the Goldreich-Julian benchmark implied 
by Eq. Q. However, since AXPs possess luminosities 
Lx 3> E SD , profound collatitudinal confinement may 
not prove necessary. Energetically, non-dipolar structure 
at the surface is easily envisaged, since the dissipation 
mechanism that powers AXP emission can restructure 
the fields; there is suggestive evidence for these reconfig- 
urations from variations seen in the pulse profiles, par- 
ticularly after flaring activity (e.g. see Kaspi et al. 2003). 
The electron energization zone may then cover a much 
larger range of colatitudes than is assigned to a standard 
polar cap, and may extend to closed field lines in equato- 
rial regions. Of course, deviations from dipole structure 
will modify the contours for the resonasphere consider- 
ably from those illustrated just below, though mostly in 
the inner magnetosphere. 



2.2 Resonant Compton Scattering 

In strong neutron star fields the cross section for Comp- 
ton scattering is resonant at the cyclotron energy and a 
series of higher harmonics (e.g. see Daugherty & Harding 
1986), effectively increasing the magnitude of the process 
over the Thomson cross section er T by as much as the or- 
der of l/(a{B) , where af is the fine structure constant. 
Here, as throughout the paper, magnetic fields are writ- 
ten in units of B cr 



i 2 e c 3 /(eh) = 4.413 x 10 13 Gauss, 
the quantum critical field strength. Klein-Nishina-like 
declines operate in supercritical fields, reducing the ef- 
fective cross section at the resonance (e.g. Gonthier et 
al. 2000). In the non-relativistic, Thomson regime (e.g. 
see Herold 1979), only the fundamental resonance is re- 
tained. For the specific case of ultra-relativistic electrons 
colliding with thermal X-rays, in the electron rest frame 
(ERF), the photons initially move mostly almost along 
B, and the cyclotron fundamental is again the only res- 
onance that contributes (Gonthier et al. 2000). 

The dominance of this resonance in forming upscat- 
tcring spectra leads to an effective kinematic coupling 
between the energies e 7 TO e c 2 and 7 e m e c 2 of colliding 
photons and electrons, respectively, and the angle of the 
initial photon # 7 to the magnetic field lines: the cy- 
clotron fundamental is sampled when 

7 e £ 7 (l — cos# 7 ) w B , for 7 e 3> 1 . (2) 

The simplicity of this coupling automatically implies that 
integration over an angular distribution of incoming pho- 
tons results in a flat-topped emission spectrum for Comp- 
ton upscattering of isotropic photons in strong magnetic 



fields. This characteristic is well-documented in the lit- 
erature in the magnetic Thomson limit (e.g. see Der- 
mer 1990; Baring 1994; for old gamma-ray burst scenar- 
ios, and Daugherty & Harding 1989; Sturner, Dermer 
& Michel 1995 for pulsar contexts), specifically for col- 
lisions between j e 3> 1 electrons and thermal X-rays 
emanating from a neutron star surface. 

It is instructive to compute the zones of influence of 
the resonant Compton process for magnetic dipole field 
geometry. For X-ray photons with momentum vector k , 
emanating from a single point at position vector R e on 
the stellar surface with colatitude 9 e (i.e., 9 e = corre- 
sponds to the magnetic pole), the kinematic criterion in 
Eq. (0) selects out a single photon angle 9 1 for a given 
local B , both of which are dependent on the altitude and 
colatitude of the point of interaction. This assumes that 
the photon propagates with no azimuthal component to 
its momentum, a specialization that will be remarked 
upon shortly. In the absence of rotation, the resonance 
criterion defines a surface that is azimuthally symmetric 
about the magnetic field axis. The locus of the projec- 
tion of this surface onto a plane intersecting the magnetic 
axis can be found through elementary geometry, assum- 
ing a flat spacetime. Observe that light bending due to 
the stellar gravitational potential will modify these loci 
significantly in inner equatorial regions. At an altitude 
r and colatitude 9 , denoted by position vector r , the 
magnetic dipole polar coordinate components are 



B r 



-Bpi? 3 



cos 9 



BpR 



(3) 



where R is the neutron star radius and B p is the surface 
field strength at the magnetic pole. At this location r , 
the corresponding polar coordinate components k r and 
kg of the photon momentum are given by 



w\ 

kg 
Ikl 



X 



V 1 - 2 XCOs( 



sin 



:) 



(4) 



^l-2 X cos(9-9 e )+ x 2 



where \ — r /R is the scaled altitude. The geometry of 
the magnetic dipole then uniquely determines the angle 
9 1 = 0-y(r/R, 6; 9 e ) of the photon to the field via the 
relation cos6> 7 = k.B/|k|.|B| = (k r B r + k g B e )/\k\.\B\ : 



cos 9~, 



2cos% - cos(fl - 9 e )} +sin6> sin(0 - 9 e 
\A + 3cos 2 6> y/l - 2xcos(<9 - 6 e ) + X 2 



.(5) 



Inserting this into Eq. (J2J) , for the specific case of 9 e = , 
yields the equation for the locus defining the surface of 
resonant scattering for outward-going electrons: 



3 _ y/l + 3 COS 2 9 
x 1 — COS &v 



Br, 



(6) 



Here <P is the key parameter that scales the altitude of 
the resonance locale, and typically might be in the range 
1 — 10 3 for magnetars of j e ~ 10 2 — 10 4 . Eq. @ can 
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Surfaces of Last Resonant Compton Scattering 




Fig. 1 Contours in a section of a pulsar magnetosphere that 
depict cross sections of the surfaces of last resonant scatter- 
ing, i.e. the maximal extent of the Compton resonasphere. 
The heavyweight contours (in dark blue) are computed for 
different values of the resonance parameter denned in 
Eq. @, and at extremely high altitudes asymptotically ap- 
proach the magnetic axis (vertical line). The filled red circle 
denotes the neutron star, whose radius R establishes the 
spatial scale for the figure. The case illustrated is for pho- 
tons emanating from the polar axis (i.e. 8 e = 0° ), denoted 
by black dots, for which the neutron star shadow regions 
are demarcated by the dotted boundaries, and only the sur- 
faces (azimuthally-symmetric about the magnetic axis) are 
accessible to resonant Compton interactions. In computing 
the contours, the upscattering electrons were assumed ultra- 
relativistic, and the depicted spatial scales are linear. 



be rearranged into polynomial form, but must be solved 
numerically. For the special case of soft photons emitted 
from the surface pole ( 9 e = ), the surfaces of resonant 
scattering for different tf' are illustrated in Fig. ^ as the 
heavy blue contours. The shadows of the emission point 
are also indicated to demarcate the propagation exclu- 
sion zone for the chosen emission colatitude. 

The altitude of resonance is clearly much lower in the 
equatorial regions, since the photons tend to travel more 
across field lines in the observer's frame, and so access 
the resonance in regions of higher field strength. In con- 
trast, at small colatitudes above the magnetic pole, 9 7 
is necessarily small, pushing the resonant surface to very 
high altitudes where the field is much lower. In this case, 
the loci asymptotically approach the polar axis at infin- 



ity, satisfying x 3 sm2 & ~ 16<? for the depicted case of 
9 e = . For 9 e ^ cases (not depicted) , the contours 
are morphologically similar, though they incur signifi- 
cant deviations from those in Fig. ^ both in equatorial 
and polar regions; for example, when 9^0, the loci 
do not extend to infinity for 9 e > . Clearly, by sam- 
pling different emission colatitudes 9 e these surfaces are 
smeared out into annular volumes. Furthermore, since 
most photons not emitted at the poles possess azimuthal 
components to their momenta, propagation out of the 
plane of the diagram must also be considered, modifying 
Eqs. JSJ and ©. In the interests of compactness, such al- 
gebra is not offered here. It suffices to observe that intro- 
ducing an azimuthal component to the photon momen- 
tum generally tends to increase propagation across the 
field in the observer's frame, i.e. increasing # 7 , so that 
the resonance criterion in Eq. J2J is realized at lower alti- 
tudes and higher field locales. Hence, taking into account 
azumithal contributions to photon propagation in the 
magnetosphere, loci like those depicted in Fig. factually 
represent the outermost extent of resonant interaction, 
and so are surfaces of last resonant scattering, i.e. the 
outer boundaries to the Compton resonasphere. It is 
evident that, for the majority of closed field lines for long 
period AXPs, this resonasphere is confined to within a 
few stellar radii of the surface. Clearly, introducing more 
complicated, non-dipolar field topologies will also tend 
to lower the altitudes of the resonasphere. 



3 Compton Upscattering Spectra in AXPs 

To gain an initial idea of what emission spectra might 
be produced in the resonasphere of AXPs, collision inte- 
gral calculations of upscattering spectra are performed. 
Here results are presented for monoenergetic electrons 
of Lorentz factor "f e , from which spectral forms for vari- 
ous electron distributions can easily be inferred. This ap- 
proach forgoes considerations of electron cooling, which 
naturally generates quasi-power-law distributions, at least 
in the magnetic Thomson limit (e.g. see Dermer 1990; 
Baring 1994); such issues will be addressed for supercrit- 
ical fields in future presentations. To simplify the formal- 
ism for the spectra, monoenergetic incident photons of 
dimcnsionlcss energy e 1 = e s will be assumed, with the 
implicit understanding that e s ~ 3kT/m e c 2 , so that val- 
ues of s s ~ 0.003 are commensurate with thermal pho- 
ton temperatures kT ~ 0.5 — 1 keV observed in AXPs 
(see Perna et al. 2001). Distributing e s via a Planck 
spectrum provides only small changes to the spectra il- 
lustrated here, serving only to smear out modest spectral 
structure at the uppermost emergent photon energies. 

Central to the characteristics of the spectral shape for 
resonant upscattering problems are the kinematics asso- 
ciated with both the Lorentz transformation from the 
observer's or laboratory frame (OF) to the electron rest 
frame (ERF), and the scattering kinematics in the ERF. 
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Since choices of photon angles in the two reference frames 
are not unique, a statement of the conventions adopted 
here is now made to remove any ambiguities. Let the 
electron velocity vector in the OF be (3 e . This will be 
parallel to B due to rampant cyclo-synchrotron cooling 
perpendicular to the field. The dimensionless pre- and 
post-scattering photon energies (i.e. scaled by m e c 2 ) in 
the OF are Si and e / , respectively, and the correspond- 
ing angles of these photons with respect to — j3 e (i.e. 
field direction) are ©i and Of , respectively. Observe 



that 04 



0~, establishes a connection to the notation 



used in Section [21 With this definition, 



cos O 



hi 



and the zero angles are chosen anti-parallel to the elec- 
tron velocity. Here, and are the initial and final 
photon three-momenta in the OF. Boosting by /3 e to 
the ERF then yields pre- and post-scattering photon en- 
ergies in the ERF of u>i and LUf , respectively, with cor- 
responding angles with respect to ~/3 e of Oi and Of . 
The relations governing this Lorentz transformation are 



'Jij = 7e£i,/(l + (3 e cosOij) 



(8) 



cos 



i,f 



COS Ojj + [j e 
1 + P e COS Oi. f 

The inverse transformation relations are obtained from 
these by the interchange Oij <-> &ij and the substitu- 
tions ujij — > Eij and e — * —fie ■ The form of Eq. JHJ 
guarantees that for most Oi , the initial scattering angle 
6i in the ERF is close to zero when j e 3> 1 , exceptions 
being cases when cos Oi — (3 e . These exceptional cases 
form a small minority of the upscattering phase space, 
and indeed a small contribution to the emergent spectra, 
and so are safely neglected in the ensuing computations. 
This Oi s=s approximation yields dramatic simplifica- 
tion of the differential cross section for resonant Comp- 
ton scattering, and motivates the particular laboratory 
frame angle convention adopted in Eq. Q. 

The scattering kinematics in the ERF differ from that 
described by the familiar Compton formula in the ab- 
sence of magnetic fields (e.g. see Herold 1979; Daugh- 
erty & Harding 1986). In the special case ftj » that is 
generally operable for the scattering scenario here, the 
kinematic formula for the final photon energy ujf in the 
ERF can be approximated by 



uj f = u'(ui, 6 S ) = 



2u>i r 



1 + \ll - 2uj t r 2 sin 2 9 f 



where 



1 + u>(l — cos I 



f) 



(9) 



(10) 



is the ratio ujf/uji that would correspond to the non- 
magnetic Compton formula, which in fact does result if 
LJiV 2 sin 2 8f <C 1. Eq. JSJ can be found in Eq. (15) of 
Gonthier et al. (2000), and is realized for the particular 



case where electrons remain in the ground state (zeroth 
Landau level) after scattering. Such a situation occurs for 
the resonant problem addressed in this paper, a feature 
that is discussed briefly below. 

Let n 7 be the number density of photons resulting 
from the resonant scattering process. For inverse Comp- 
ton scattering, an expression for the spectrum of pho- 
ton production dn y /(dt de f dfif) , differential in the pho- 
ton's post-scattering laboratory frame quantities e f and 
li f = cos Of , was presented in Eqs. (A7)-(A9) of Ho and 
Epstein (1989), valid for general scattering scenarios, in- 
cluding Klein-Nishina regimes. This was used by Dermer 
(1990) and Baring (1994) in the magnetic Thomson do- 
main, i.e. when B < 1 and the photon energy in the 
electron rest frame (ERF) is far inferior to m e c 2 . Such 
specialization is readily extended to the magnetar regime 
by incorporating into the Ho and Epstein formalism the 
magnetic kinematics and the QED cross section for fully 
relativistic cases. The result can be integrated over fif 
and then written as 



da 



(11) 



dtdfj-H^J^ 

S[ Uf - U '(u u 9 t )] 7e (^ e %) die^ej) • 

noting that the angle convention specified in Eq. J7J re- 
quires the substitution (3 e — > — f3 e in Eqs. (A7-A9) of Ho 
and Epstein (1989). Here, the notation \ii — cos (9,: and 
[if — cos Of is used for compactness, n e is the num- 
ber density of relativistic electrons, and n s is that for 
the soft photons. These incident monoenergetic photons 
are assumed to possess a uniform distribution of angle 
cosines fii in some range [i^ < /Ltj < /i+ , which is gen- 
erally broad enough to encompass the resonance, i.e. the 
value /li = [B/{~f e e s ) — l]//? e ■ The bounds on the (if in- 
tegration, defining the observable range of angle cosines 
with respect to the field direction, will be specialized to 
fil = — 1 and fi u = 1 in the illustration here, though 
dependence of the emergent spectra Of values will be 
discussed below. The function w'(wi, Of) , which appears 
in the delta function in Eq. (|ll|l that encapsulates the 
scattering kinematics, is that defined in Eq. 

Fully relativistic, quantum cross section formalism 
for the Compton interaction in magnetic fields can be 
found in Herold (1979), Daugherty & Harding (1986), 
and Bussard, Alexander & Meszaros (1986). These ex- 
tend earlier non-relativistic quantum mechanical formu- 
lations such as in Canuto, Lodenquai & Rudcrman (1971), 
and Blandford & Scharlemann (1976). The differential 
cross section, do jd cos Of , appearing in Eq. (|ll|l is taken 
from Eq. (23) of Gonthier et al. (2000), and incorpo- 
rates the relativistic QED physics. Yet, it is specialized 
to the case of scatterings that leave the electron in the 
ground state, the zeroth Landau level that it originates 
from. This expedient choice is entirely appropriate, since 
it yields the dominant contribution to the cross section 
at and below the cyclotron resonance (e.g. Daugherty & 
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Harding 1986; Gonthier et al. 2000). For considerations 
below, where information on the final polarization state 
of the photon is retained, Eq. (22) of Gonthier et al. 
(2000) is used; this simplifies to the forms 



da 



d cos f 



\,X 3a T 0"/) 3 T\ U exp{-u;}sm 2 e f /[2B}} 



32 



[1 + u>i(l - cos9f) - LUf sin Of] 



1 



1 



'12) 



{uji - By + (r cyc /2y (ui + B-cy 



(13) 



for the differential cross sections, where the cyclotron de- 
cay width r cyc (discussed below) has been introduced to 
render the resonance finite in the form of a Lorentz pro- 
file. The exponential factor in Eq. Ijl2(l is a relic of the La- 
guerre functions that signal the discretization of momen- 
tum/energy states perpendicular to the field. Further- 
more, the factor in square brackets in the denominator 
is always positive, being proportional to l — Uf-r sin 2 Of , 
which can be shown to be precisely the square root ap- 
pearing in Eq. ©. Also, £ = u;,o;/(l — cosOf) , and 

Tj| = 2 cos 2 Of +Wj(l -cos6>/) 2 -Lv f sm 2 f 

T± = 2 + cj 4 (1 - cosOf) 2 -ujf sin 2 Of . 

Here, the standard convention for the labelling of the 
photon linear polarizations is adopted: || refers to the 
state with the photon's electric field vector parallel to 
the plane containing the magnetic field and the photon's 
momentum vector, while _L denotes the photon's electric 
field vector being normal to this plane. 

For magnetic Compton scattering, in the particular 
case of photons propagating along B prior to scattering 
(i.e. 0i = 0), the differential cross sections are indepen- 
dent of the initial polarization of the photon (Gonthier 
et al. 2000); this property is a consequence of circular po- 
larizations forming the natural basis states for 0i = . 
Therefore, transitions _l_ — »- 1 1 and ||— >|| yield identical 
forms for the cross sections, and separately so do the 
transitions _L— >J_ and ||— >-L . Accordingly, the cross sec- 
tions in Eq. I|12l) are labelled only by the post-scattering 
linear polarization state of the photon; they are summed 
when polarization-independent results are desired, i.e. 
da/dcosOf — (da\\/dcos0f + da±/dcos0f) . Therefore, 
clearly the upscattered photon spectra presented here 
are insensitive to the initial polarization level (zero or 
otherwise) of the soft photons. 

The integrals in Eq. i|ll[) can be manipulated using 
the Jacobian identity d^d^if — du>i duif / ( r y 2 f3 2 Sie /) to 
change variables to ioi and u)j . Observe that the values 
of fi± do not impact these integrations provided that 
the resonance condition in Eq. @ is sampled. The ujf 
integration is then trivial. The Ui integration is more 
involved, but can be developed by suitable approxima- 
tion, as follows. The relativistic Compton cross section is 
strongly peaked at the cyclotron fundamental (see Fig. 2 
of Gonthier et al. 2000) due to the appearance of the 
resonant denominator l/[(w, — B) 2 + (r cyc /2) 2 ] , where 
r cyc <C B is the dimensionless cyclotron decay rate from 



the first Landau level. Therefore, this Lorentz profile can 
be approximated by a delta function in space of iden- 
tical normalization in an integration over u>i : 



1 



(w< - Bf + {r cyc /2y 



2tt 



eye 



S(ui - B) 



(14) 



This mapping, which was adopted in the magnetic Thom- 
son limit by Dcrmer (1990), renders the uJi trivial, with 
the non-resonant term in the cross section in Eq. (|12|) 
being neglected, and the evaluation of the integrals in 
Eq. (jl 1|> complete. The spectra scale as the inverse of the 
decay rate r cyc , whose form can be found, for example, 
in Eqs. (13) or (23) of Baring, Gonthier & Harding (2005; 
see also Latal 1986; Harding & Lai 2006). For B <C 1 , 
r cyc w 4a{B 2 /3 , which traces classical cyclotron cooling, 
while for B ^> 1 , quant um ef fects and recoil reductions 
generate r cyc w (ctt/e)y/B/2. 
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Fig. 2 Resonant Compton upscattering spectra (scaled) 
such as might be sampled in the magnetosphere of an AXP, 
for different relativistic electron Lorentz factors 7 e , as la- 
belled. The emergent photon energy Sf is scaled in terms 
of m e c 2 . The chosen magnetic field strengths of B = 3B CV 
(heavyweight, blue) and B = 0.3B C r (lighter weight, red) 
correspond to different altitudes and perhaps colatitudes. Re- 
sults are depicted for seed photons of energy e s ~ 0.003 
(marked by the green vertical line), typical of thermal X- 
rays emanating from AXP surfaces; downscattering resonant 
emission at e/ < e s was not exhibited. 



Representative spectral forms are depicted in Fig. [3 
for the situation where the emergent polarization is not 
observed. Because of the approximation to the resonance 
in Eq. (|14|) . non-resonant scattering contributions were 
omitted when generating the curves; these contributions 
produce steep wings to the spectra at the uppermost and 
lowermost energies (not shown), and a slight bolstering 
of the resonant portion. This resonant restriction suffices 
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for the purposes of this paper, and kincmatically limits 
the range of emergent photon energies et to 

7e (l + /3 e )ff 

1 + 2B 



7e(l-/3 e )B < Ef 



< 



(15) 



This range generally extends below the thermal pho- 
ton seed energy e s . Such downscatterings correspond 
to forward scatterings cos Of « 1 in the ERF such that 
Of 6i ~ l/7e • Since they constitute a miniscule por- 
tion of the angular phase space (and energy budget), 
only upscattering spectra are exhibited in the Figure. 

In the B = 0.3 case, a quasi- Thomson regime, the 
spectra show a characteristic flat distribution that is in- 
dicative of the kinematic sampling of the resonance in 
the integrations (e.g. see Dermer 1990; Baring 1994). For 
much of this spectral range, forward scattering in the 
observer's frame is operating: /if = cos (9/ « 1 . This es- 
tablishes almost Thomson kinematics, with the scattered 
photon energy in the ERF satisfying u>f w u>i . Substi- 
tuting these approximations into Eqs. (fT2*)l and lfT3"|l. and 
then summing over final polarizations yields a total ap- 
proximate form for the flat portions of the spectral pro- 
duction rate in Eq. (|llf> for cases 7 e 1 : 



n e n s a T c 



3ttB 2 



dtdef ~ 4T cyc7e ^ ' (16) 

The magnetic field and 7 e dependences are evident in 
Fig. |21 (remembering that the curves are multiplied by 
7g for the purposes of illustration), and since r cyc oc B 2 
when B <C 1 , the normalization of the flat spectrum is 
independent of the field strength in the magnetic Thom- 
son regime. Only at the highest energies does the spec- 
trum begin to deviate from flat (i.e. horizontal) behav- 
ior, and this domain corresponds to significant scattering 
angles in the ERF, i.e. cosines 1 — cos Of not much less 
than unity. Then the mathematical form of the differen- 
tial cross section becomes influential in determining the 
spectral shape. Specifically, for Of ~ it/2, Tu drops far 
below T_l , as is evident in Eq. I|13fl , causing the observed 
dip in the spectra. At slightly higher energies e / , there 
is a recovery when Tj| rises as Of — > tt . Observe that Tj_ 
is far less sensitive to scattering angles Of in the ERF, 
except for supercritical fields when recoil becomes signif- 
icant. The sum of the two contributions yields a slight 
cusp at the maximum energy Ef « 7 e (l + (3 e )B/(l + 2B) 
when B ^ 0.5 , which disappears for higher fields when 
the recoil reductions of da/d cos Of become dominant. 

In AXPs, the B = 0.3 case best represents higher 
altitude locales for the resonasphere, such as at smaller 
colatitudes near the polar axis. For B = 3 , more typi- 
cal of equatorial resonance locales, the flat spectrum still 
appears at energies e s < £/ <C 7e(l + @e) B/(l + 2B) , 
when again cos Of 1. Yet the curves in Fig. [21 dis- 
play more prominent reductions at the uppermost ener- 
gies Ef ~ 7e(l + (3 e )B/{l + 2B) due to the sampling 
of 1/2 < ujf <C Wj values in the ERF that correspond 
to strong electron recoil effects. Photons emitted in this 
regime have 1 — cos Of ~ 1 in the ERF and are highly 



beamed along the field in the observer's frame, as will be- 
come apparent shortly. At these maximum energies, the 
approximations cos Of » 1 and u>f « B/(l + 2B) yield 
an analytic result for the emission rate, precisely Eq. (|16|) 
multiplied by the factor 1/(1 + 2B) 2 that controls the 
severity of the reduction at the uppermost resonant en- 
ergies. Note that Klein-Nishina reductions in the cross 
section regulate the overall normalizations of the B = 3 
curves, as does the fact that the cyclotron width r cyc 
no longer scales with field strength as B 2 . 

The dimensionless electron energy loss rate d"f e /dt 
can be obtained by multiplying the differential spectrum 
dn^/dtdef in Eq. 1)11(1 by Et/n e and integrating over 
e f . Because of the flat nature of the spectrum, this re- 
ceives a dominant contribution from near the maximum 
upscattered energy 7 e (l + (3 e )B/(l + 2B) . Then, us- 
ing Eq. (|16|l . in the magnetic Thomson limit B <C 1 , 
since energies £/ ~ 2j e B contribute most, the result 
dj e /dt oc B^/r cyc oc B 2 is obviously realized for r cyc « 
4afi? 2 /3. The full integration of Eq. (fTTJ) over Ef (or 
equivalently cos Of ) is analytically tractable, leading to 



die 

dt 



n s a T c 3irB 



B < 1 



(17) 



This result is commensurate with the form derived in Eq. 
(24) of Dermer (1990). In contrast, when B 1 , the 
overall normalization of the spectrum at energies £/ <S 
7 e (l + P e )B/(l + 2B) is still controlled by Eq. JT^, but 



now the cyclotron decay rate dependence r c . 



oc 



£1/2 



and recoil reductions at the highest ej come into play, 
so that the cooling rate d-f e /dt possesses a dependence 
on the field strength that is much weaker than £? 3 / 2 . 

The resonant upscattering spectra are potentially po- 
larized, perhaps strongly. Isolating the specific polariza- 
tion forms in Eq. (|13fl . the polarization-dependent reso- 
nant Compton spectra are readily computed and are il- 
lustrated in Fig. |21 It is clear from Eq. (|13|) that the emis- 
sivities for photons of final polarization state _L should 
always be superior to those for the || state. In a clas- 
sical description, this is a consequence of the physical 
ease with which an oscillating electron can resonantly 
drive emission with electric field vectors perpendicular 
to B. Yet for much of the range of emergent photon en- 
ergies above e s — 0.003 , there is no material difference 
between fluxes for the two final polarizations. This case 
corresponds to cos Of w 1, i.e. photon emission along 
the field in the ERF, and hence induces zero linear po- 
larization by symmetry, but permits emergent circular 
polarization. At the highest produced energies, signifi- 
cant differences between _L and || emission appear, with 
1 — cos Of no longer very small (see Eq. i|13|) to help 
identify this characteristic). This domain motivates the 
development of medium energy gamma-ray polarimeters 
as a tool for geometry diagnostics. 

Another feature of the upscattering process that is 
highlighted in the Fig. |31 is the intense beaming of ra- 
diation along the field in the observer's frame of refer- 
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Fig. 3 Resonant Compton upscattering spectra appropri- 
ate for an AXP magnetosphere, scaled as in Fig. |^1 but this 
time displaying the two polarizations _L (heavyweight) and 
|| (lightweight) for the produced photons; the dominance of 
the J_ polarization near the uppermost energies is evident. 
Again, the emergent photon energy £/ is scaled in terms 
of m e c 2 , and results are presented for different relativis- 
tic electron Lorentz factors 7 e , as labelled. Specific emer- 
gent angles of the emission in the observer's frame, with re- 
spect to the magnetic field direction, are indicated by the 
filled magenta symbols, with triangles denoting Of — 5° , 
squares corresponding to Of = 1° , and circles represent- 
ing Of — 0.2° . The magnetic field strengths of B = 3B CY 
(blue) and B — 0.3B CI (red) correspond perhaps to lower 
and higher altitudes, respectively. Again, the green vertical 
line marks the soft photon energy e 3 = 0.003 , with down- 
scattering resonant emission at £/ < e s not being exhibited. 

ence, and the profound correlation of the angle of emis- 
sion O f with the emergent photon energy e / . This tight 
correspondence is mathematically guaranteed by the ap- 
pearance of the kinematic delta function in Eq. (|llfl to- 
gether with the intrinsically narrow nature of the scat- 
tering resonance, which permits the delta function ap- 
proximation in Eq. I)14|l . The extremely narrow range of 
j for each observed energy £/ is broadened when non- 
resonant scattering is introduced. In the resonant case 
here, the 7 e 3> 1 regime dictates that most of the emis- 
sion is collimated to within 5° of the field direction, 
and rapidly becomes beamed to within 0.2° as the fi- 
nal photon energy increases towards its maximum. This 
kinematic characteristic guarantees that spectral forma- 
tion in Compton upscattering models is extremely sen- 
sitive to the observer's viewing perspective in relation 
to the magnetospheric geometry, offering useful probes 
if pulse-phase spectroscopy is achievable. 

It should be noted that the spectra in Figs. El and 
are in principal subject to attenuation by magnetic pair 
production, 7 — > e + e~ , reprocessing the highest energy 
photons to lower energies. This process is sensitive to the 



angle Of of the scattered photons to the field, which, 
from the results depicted here, is strongly coupled to 
their energy Sf . Gonthier et al. (2000) demonstrated 
(see their Figure 7) that generally, for an extended range 
of values for B , pair creation attenuation would only 
operate for uji <; 10-30 in the ERF, since the resonant 
Compton process couples the values of ujf and 6f . For 
the resonant scattering considerations here, this crite- 
rion translates to 7 — + e + e~ being rife for local fields 
B i; 10 and being marginal, or more probably ineffec- 
tive, at lower field strengths. For B ^ 10 scattering 
circumstances, any pair creation that ensues does so by 
generating pairs in the ground state (e.g. Usov & Mel- 
rose 1995; Baring & Harding 2001), with Lorentz factors 
less than j e since ej < j e . Accordingly the cascading 
would consist at first of generations of upscattering and 
subsequent pair creation, until high enough altitudes are 
encountered for pair production to access excited Lan- 
dau levels for the pairs, and then synchrotron/cyclotron 
radiation can ensue and complicate the cascade. 



4 Discussion 

It is clear that the spectra exhibited in Fig. [21 are con- 
siderably flatter than the hard X-ray tails (<~ el 1 ) seen 
in the AXPs, and extend to energies much higher than 
can be permitted by the Comptel upper bounds to these 
sources. However, they represent a preliminary indica- 
tion of how flat the resonant scattering process can ren- 
der the emergent spectrum, and what an observer de- 
tects will depend critically on his/her orientation and 
the magnetospheric locale of the scattering. The one-to- 
one kinematic correspondence between Sf and fif (il- 
lustrated via the filled symbols in Fig. imposed by 
ujf = 7 e £/(l — fif) = w'(aii) with LU{ = B , implies that 
the highest energy photons are beamed strongly along 
the local field direction. This may or may not be sampled 
by an instantaneous observation, which varies with the 
rotational phase. Realistically, depending on the pulse 
phase, angles corresponding to /1 f < 1 will be predomi- 
nant, lowering the value of £/ . Yet how low is presently 
unclear, and remains to be explored via a model with full 
magnetospheric geometry, an essential step. One can also 
expect substantial spectral differences between scatter- 
ing locales attached to open and closed field lines, and 
also between dipolar and more complicated field mor- 
phologies with smaller radii of curvature. 

Distributing the electron j e such as through reso- 
nant cooling will generate a convolution of the spec- 
tra depicted in Fig. |2 observe that the 7^ scaling of 
the y-axis implies that the normalization of the curves 
is a strongly-declining function of 7 e . This can clearly 
steepen the continuum for a particular range of fif . For 
example, since £/ oc 7 e near the maximum photon en- 
ergy for a fixed electron Lorentz factor, integration over 
a truncated r y~ p distribution, where 7 e > 7 e .mm , natu- 
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rally yields a photon spectrum e ^ at energies above 
the critical value Sf ~ 2j e ,minB / (1+25) where the res- 
onant flat top turns over. In particular, cooling the elec- 
trons as they propagate in the magnetosphere can lead 
to significant and possibly dominant contributions from 
Lorentz factors j e ^5 10 at higher altitudes and lower B 
that can evade the Comptel constraints on the AXPs (see 
Kuiper et al. 2006). In the Thomson regime, the cooling 
tends to steepen the continuum (e.g. Baring 1994) in 
the X-ray band due to a pile-up of electrons at low j e . 
If electrons propagate to high altitudes in magnetars, 
a similar steepening should be expected. As the polar- 
ized signal appears at the highest energies for each -f e , 
a somewhat broad range of energies will exhibit polar- 
ization above around 50-100 keV when integrating over 
an entire cooled electron distribution. Note also that the 
cooling may persist down to mildly-relativistic energies, 
i.e., 7 e ~ 1 , in which case it can seed a multiple scatter- 
ing Comptonization of thermal X-rays that may gener- 
ate the steep non-thermal continuum observed in AXPs 
below 10 keV. Such resonant, magnetic Comptonization 
has been explored by Lyutikov & Gavriil (2006), and 
provides very good fits to both Chandra (Lyutikov & 
Gavriil 2006) and XMM (see Rea et al. 2006) spectral 
data for the AXP IE 1048.1-5937. 

Finally, the introduction of non-resonant contribu- 
tions can have a significant impact on the spectral shape. 
The resonasphere is spatially confined, almost to a sur- 
face for a given photon trajectory, so that for most of an 
X-ray photon's passage from the stellar surface, it scat- 
ters only out of the cyclotron resonance with outward- 
going electrons, particularly if rapid cooling is operat- 
ing (a common circumstance). Moreover, inward-moving 
electrons traversing closed field lines participate in head- 
on collisions with surface X-rays, and so have great diffi- 
culty accessing the resonance. These circumstances pro- 
vide ample opportunity for the emergent spectrum to 
develop a significant non-resonant component that does 
not acquire the characteristically flat spectral profiles 
exhibited here. Since Eq. I|15|) establishes £/ 2j e B 
for B <C 1 , then this upper bound becomes inferior 
to the classical, non-magnetic inverse Compton result 
Ef ~ 47g£ s /3 when B -f e e s . This then defines a global 
criterion for when non-resonant Compton cooling domi- 
nates the resonant process. Assessing the relative weight 
of the resonant and non-resonant contributions requires 
a detailed model of magnetospheric photon and electron 
propagation and Compton scattering; this will form the 
focal point of our upcoming modeling of the high energy 
emission tails from Anomalous X-ray Pulsars. 
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